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I.  INTRODUCTION 


The  accuracy  of  geoid  heights  and  vertical  deflections  derived  from  satellite  altimetry  depends  pri- 
marily on  the  accuracy  of  the  satellite’s  ephemeris  and  on  the  measurement  noise.  The  measurement  noise 
is  for  most  part  a .^lort-wavelength  (<  1 km)  random  process;  it  is  usually  taken  out  by  filtering  or  smooth- 
ing the  along-track  geoid  heights.’  The  ephemeris  errors^  on  the  other  hand,  consist  mainly  of  long- 
wavelength  (>  1000  km)  biases  and  bias  rates.  These  can  be  solved  for  and  eliminated  if  a sufficient  number 
of  track  intersections  is  available  in  a given  region.  The  data  for  this  solution  are  the  geoid-height  discrep- 
ancies (differences  of  the  two  filtered  geoid  heights)  at  the  intersections. 

This  report  presents  an  approximate  least-squares  algorithm  for  such  a solution.  It  is  an  iterative 
algorithm  in  the  sense  that  first  it  solves  for  biases  alone.  Subtracting  the  effect  of  this  solution  from  the 
original  (Oth-order)  discrepancies  leaves  a set  of  Ist-order  discrepancies.  These  are  then  used  to  solve  for  bias 
rates,  leaving  2nd-order  discrepancies,  which  are  used  to  solve  for  2nd-order  effects,  etc.  This  process  can  be 
terminated  at  the  user’s  discretion. 

It  is  easy  to  show  that  this  iterative  procedure  does  not  give  the  optimum  least-squares  solution.  For 
example,  if  the  algorithm  is  terminated  with  the  second-order  solution,  the  result  will  not  be  as  good  as  a 
least-squares  solution  with  the  full  parameter  set  consisting  of  biases,  bias  rates  and  2nd-order  terms.  If 
biases  alone  are  desired,  then  this  algorithm  of  course  will  give  the  correct  least-squares  solution.  It  is  dif- 
ficult to  quantify  by  how  much  the  iterative,  segmented  solution  is  degraded.  Test  cases  show  that  for  most 
applications  this  will  be  negligible.  It  has  the  great  advantages,  however,  that  it  is  relatively  simple  and  re- 
quires much  less  computer  time  and  storage.  This  comes  from  the  fact  that  the  matrices  to  be  inverted  in 
this  algorithm  have  dimensions  which  depend  only  on  the  number  of  tracks,  and  not  on  the  number  of 
tracks  times  the  number  of  bias  parameters  to  be  solved  for,  as  would  be  the  case  with  the  full  least-squares 
solution. 

The  general  algorithm,  which  is  derived  in  chapter  2,  takes  account  of  missing  data  points  (where  no 
intersections  exist)  by  giving  them  zero  weight.  All  other  data  are  weighted  equally.  The  problem  is 
inherently  singular  with  infinitely  many  solutions.  The  singularities  are  lifted  by  assigning  a-priori  variances 
for  the  bias  parameters  on  each  track.  The  track  biases  ate  assumed  to  be  uncorrelated.  In  chapter  3 the 
solution  is  presented  for  biases  only  when  there  are  no  missing  data  points.  This  solution  turns  out  to  be 
very  simple;  no  matrix  inversions  are  needed. 

In  chapter  4,  several  test  cases  are  presented  which  will  be  useful  for  check-out  purposes.  Finally, 
Appendix  A gives  a step-by-step  procedure  for  coding  the  algorithm  on  a computer. 


II.  DERIVATION  OF  THE  ALGORITHM 
A.  THE  observation  EQUATIONS 

Figure  1 shows  a network  of  w + m intersecting  tracks  which  consist  of  n “rows,”  1 < / < n.  and  m 
“columns,”  I <i<m.  The  index  i will  always  denote  a row,  and  / a column. 

We  thus  have  a total  of  n'm  intersections.  At  each  of  these  intersection  points  (i,  /)  we  cither  have  a 
data  point  Ajf,  or  the  data  point  is  missing  in  which  case  an  asterisk  is  placed  at  that  intersection.  (In 
practice  a missing  data  point  means  that  there  is  no  intersection.) 


1 


Figure  1 . Intersection  Data  A,y 


The  model  for  the  biases  on  each  track  (row  or  column)  is  assumed  to  be  a power  series  in  t - the 
time  along  the  track,  with  r = 0 at  the  middle  of  the  track: 

6(t)  = ft  + ftr  + + ... />***T*  (1 

We  shall  first  solve  for  the  constant  biases  alone.  Subtracting  their  contribution  from  the  original  data 
will  give  a new  set  of  residuals  MjJK  This  new  set  of  data  will  then  be  used  to  solve  for  the  bias  rates  ft, 
yielding  2nd-order  residuals  etc.,  until  any  desired  order  k. 

The  observation  equations  for  a least-squares  solution  of  any  order  k are  given  by 

A*.*'  =x  lJ  ■ - y V- 
^1/  '*1*^/1  -I  II' 

1 < /■  < w,  (2 

I </<  w. 


where  the  unknowns  x,-,Vy,  and  the  constants  U/j,  Vjj  are  defined  in  Table  1 . 


The  times  r in  this  table  are; 


= jih  tifng  on  /*  row, 
T,-^^  ~ tin’®  on  / * column. 


Table  1.  The  Unknowns  and  Ujj,  V^j  of  Equation  (2) 


Solution  Order 

^0- 

0 (bias) 

1 

1 

1 (bias  rate) 

rW 

/' 

2 

Irfl^ 

k 

ifY 

B.  THE  NORMAL  EQUATIONS 

The  normal  equations  for  the  least-squares  solution  of  Equation  (2)  with  arbitrary  weighting  W are 


A^WA  (^'^  = A'^WA, 


(3) 


where  A is  the  (n’m)X  1 data  vector  formed  from  the  residuals  A(^l.  A is  the  (n'm)X(n  + m)  matrix  of 
partial  derivatives 

3A* 


dy,- 


(4) 


Kk^n-m,  l<i<w,  l</<»i. 


It  is  written  out  explicitly  in  Table  2.  The  weighting  matrix  W will  be  used  here  only  to  accommodate 
missing  data  points,  although  it  can  be  implemented  in  general  as  an  arbitrary  weighting  matrix  for  all  of 
the  data.  All  existing  data  Ajf  will  be  given  unit  weight.  Accordingly,  let  W be  the  (n"wi)X(n*m)  diagonal 
matrix  with 


= 1 fork  = (/, /)  where  A,-,-  exists, 

= 0 for  k = (t.  /)  where  Ajy  is  missing. 
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Table  2.  Theyl  Matrix 


0 0 (/,3 
0 0 f/23 

0 0 f/33 


-'"2m 


It  is  easy  to  show  that  this  leads  to 

I 

i A^h'  = A^ 

i 

^ providing  that  we  redefine  f/y,  and  K,y  such  that 


(5) 

(6) 


at  all  intersections  (i,  j)  with  missing  data  points,  and  leave  them  as  given  in  Table  I for  points  (/,  /)  where 
Aji  does  exist.  Equation  (3)  now  becomes 


C.  THE  SOLUTION 


The  matrix  A in  Equation  (7)  is  singular.  For  the  bias  solution  (^  = 0 in  Equation  (2)),  the  reason 
for  that  singularity  is  obvious:  we  may  add  to  ail  of  the  biasesx,  and  ty  in  Equation  (2)  an  arbitrary 
constant  without  changing  the  data  vector  Ay.  For  the  higher-order  solutions  similar  singularities  can  be 
established.  The  problem,  therefore,  has  infinitely  many  solutions.  In  order  to  obtain  a unique  solution 
we  lift  these  singularities  by  imposing  a-priori  constraints  a,-  and  (•  for  the  row  and  column  biases, 
respectively.  The  solution  of  Equation  (7)  then  is 


= B-^A^A 


with 


B = AT'A  + S, 

where  2 is  the  diagonal  matrix  of  the  reciprocal  a-priori  variances  with  elements 


kl  „2 


hr 


(8) 


(9) 


(10) 


With  the  A matrix  as  given  in  Table  2,  it  is  easy  to  obtain  the  B matrix  as 


B = 


I ^1 

I 
I 


j ^WXw) 


T ' 

^liiiXn)  I ^hmXm) 


where  A and  M are  diagonal  matrices 


Mj, 


and 


°f.R 


(II) 


(12a) 


(12b) 


(12c) 


The  range  of  indices  in  Equations  ( 1 2)  are  as  shown  by  the  dimensions  of  the  four  matrices  in  Equation  (II). 
The  inverse  of  the  B matrix  is 

'^W  X\ 


S''  = 


y Z 


(l.H 


S 


W = {^-Q^^^Q^')-\  (13a) 

Z = (Af (13b) 

X=-\-^QZ,  (13c) 

Y=-^^^Q'^W.  (13d) 

Note  that  ) = , but  we  prefer  to  leave  them  in  this  form  for  symmetry  reasons.  The  dimension  of  the 

largest  matrix  to  be  inverted  (Iv'  or  7.)  is  either  (/iX/i)  or  (wXw),  the  number  of  rows  or  columns.  These 
inverses  will  in  general  have  to  be  done  on  the  computer.  In  the  case  there  are  no  missing  data  points,  they 
can  be  done  analytically  for  the  Oth-order  (biases)  solution  (see  the  next  chapter). 

We  may  write  the  solution.  Equation  (8),  in  a little  more  expanded  form.  The  (wX  I ) Row  Solutions 


and  the  (wX  1 ) Column  Solutions  are 


X = 


V = 


= A^A. 


i fi  \ i 

The  vectors  s and  s'  are  appropriately  weighted  sums  of  the  data  residuals  on  each  row  and  column. 


s)«>=  £ !</<„; 

/=! 

n 


Equations  ( 14)  with  ( 1 2),  ( 1 3)  and  ( 1 6)  constitute  the  complete  solution  in  matrix  form.  In  Appendix 
A the  algorithm  is  given  in  component  form,  which  is  suitable  for  direct  coding  on  a computer. 


III.  SOLUTION  FOR  BIASES  WITH  A FULL  DATA  SET 

Wlien  there  arc  no  missing  data  points  Ajj  the  solution  for  biases  alone  can  be  obtained  in  closed  form, 
i.e.,  the  B matrix.  Equation  (II),  can  be  inverted  in  closed  form.  This  is  a big  advantage  over  the  general 
solution,  so  that  in  applications  where  only  the  constant  biases  are  needed,  it  may  be  advisable  to  Till  in  the 
missing  data  points. 


With  a full  data  set,  for  the  bias  solution. 


4 


(17) 


so  that  Equation  ( 1 2c)  becomes 


(18) 


which  may  be  written  as 


Q = -uv‘. 


(19) 


where 


l'^ 

1 

1 

, V = 

vl 

XI 

(20) 


are  the  (nX  1)  and  (mX  1)  constant  vectors  whose  components  are  all  unity.  We  have  the  relations 


u7’u  = n,  = 


(21) 


We  specialize  this  derivation  a little  as  compared  with  the  general  solution,  by  not  allowing  arbitrary 
a-priori  variances.  Rather,  we  require  that  all  rows  have  the  same  a-priori  variance  and  all  columns 
another,  a^.  The  more  general  case,  with  different  a’s  for  each  track,  can  be  handled  in  closed  form  also, 
but  requires  much  more  algebra.  Equations  ( 1 2a)  and  ( 1 2b)  now  become 


A = X/, 


nXn’ 


mXm' 


(22) 


X = m — , 11  = n + — ; 


so  that  Equations  (13a)  and  (13b)  can  be  written  as 


(23a) 


(23b) 
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These  matrices  can  be  inverted  by  using  the  general  identity  for  any  two  vectors  a and  b 


(/>ab^)-'  =/  + 


ab^ 

1 ± a^b 


The  result  is 


Kquations  ( 1 3c)  and  ( 1 3d)  yield 


1^=— (/+  — uu^ 

XV  a 


Z=-f/  + -^vv7’), 

d 


(f  = Xfj  - mn. 


A'=  — uv^;  K=  — vu^. 
d d 


The  solution.  Equations  (14),  can  now  be  readily  evaluated  to  be 
Row  Biases: 


/j(^)  = _L  n(R)  -I-  ^ e J.  — L C 1 < I < II 


Column  Biases: 


where  Equations  ( 16)  give 


and  where  we  have  defined 


1 (0  + 5 +-L^  l</< 

' I nd  ^ d ^ 


i=\  /=i 


■v«E>r  -v-E*;''’ 

i=  1 /=  1 


m 


(24) 


(25a) 

(25b) 

(25c) 

(26) 


(27a) 


(27b) 


(28) 


(2d) 


i 

I 


Note  that  the  full  data  set  A,y  must  satisfy  the  constraint 


S = S,f  =-S(-. 


With  this,  the  solution.  Equation  (27),  assumes  a particularly  simple  form: 


Row  Biases: 


Column  Biases: 


= -i  M«)  . A_\  , 1 

' ^ V' 


*p  = _LL(C)+, 


1 </ < w. 


IV.  NUMERICAL  EXAMPLES 

We  present  here  several  numerical  examples.  These  may  be  used  as  test  cases  for  checking  the 
algorithm  out  after  it  is  programmed  on  a computer.  The  first  four  have  only  six  intersections,  so  that  the> 
are  simple  enougli  for  hand  calculations.  The  solutions  as  well  as  some  intermediate  restdts  are  given.  The 
last  two  examples  consist  of  50  rows  and  60  columns. 


A.  EXAMPLE  1 

« = 3 rows,  m = 2 columns. 

A-priori  statistics:  = 3. 

We  construct  data  for  this  example  by  assuming  that  the  five  tracks  have  the  following  true  biases: 


°TRUi:  \ ' TRUE  to 


The  data  matrix  = h.  - bj  therefore  is 


A<'»=  -7  -2  , 
-4  I 


which  has  a standard  deviation 


o,f„  = 4.535. 


<■) 


Sl'me  intermediate  results  (see  Equations  (22)  and  (30))  are 

\ = H = S = .5.  (35) 

The  solution,  Equation  (31 ),  is: 


which  does  not  compare  well  at  all  with  the  true  biases  given  in  Equation  (32).  The  solution,  of  course,  is 
completely  insensitive  to  an  overall  bias,  simply  because  the  data  A®'  are  so  also.  We  must  compare,  there- 
fore, Equation  (36)  with  the  true  biases  minus  their  mean  (which  isb  = 2): 

“'Suu-S-^) 

Doing  that,  we  find  that  the  solution  errors  are  -5.5,  -5.0,  -4.2,  -3.3  and  -3.9%,  respectively  for  the  five 
biases.  This  implies  that  our  chosen  a-priori  a was  too  small. 

The  Ist-order  residuals,  formed  by  subtracting  the  solution.  Equation  (36)  from  the  data.  Equation 
(33),  are: 

.121  .299 

,^(1)=  . 300  -.122  , (38) 

-.142  .036 

which  has  the  standard  deviation 

V')  = 0“'4- 

This,  when  compared  with  Equation  (34),  shows  the  relative  efficacy  of  the  solution  to  reduce  the  inter- 
section discrepancies. 

B.  EXAMPLE  2 

This  is  the  same  as  Example  1 , except  that  the  a-priori  variances  arc  increased  to 

— 1 0. 

'0 


j 


I 


The  solution  is 

/ 3.979  \ 
=[  -3.981 
^-.996  / 

which  now  agrees  much  better  with  Equation  (37). 
The  remaining  Ist-order  residuals  are 


b<C) 


= { 2-991  \ 
1.993^  ’ 


(40) 


a“>  = 


.012  .028 
.028  -.012 


.013 


.003 


with  the  standard  deviation 


C.  EXAMPLE  3 


a^(,)  = .020. 


(41) 


(42) 


Here  we  take  the  same  data  as  in  the  two  previous  examples,  but  delete  the  intersection  of  the 
second  row  with  the  first  column.  We  keep  the  a-priori  variances  as  in  Example  2: 

Because  there  is  now  a “missing  data  point,”  the  general  solution  of  section  IIC  must  be  used. 

The  data  matrix  is 


A(0)  = 


1 6 
♦ _2 
-4  I 


(43) 


where  the  asterisk  indicates  the  missing  data  point;  it  can  be  filled  in  with  an  arbitrary  value.  Equations  (4) 
and  (12)  become 


(44a) 


u= 

1 0 r 

1 I 1 

, y= 

1 r 
0 1 

j,  Q = - 

1 1 

0 1 

1 

r_  - . 

_ “1 

.1  1. 

. 1 1. 

A = 


2.01  0 


0 


0 1.01  0 


_ U • '-I 

r~ 

, M = 

2.01 

0 

0 

3.01 

L 

(44b) 


Equations  (1 3)  yield 


(vrl  = 


1.180262  -.332226  -.829738 
-.332226  .677774  -.332226 

-.829738  -.332226  1.180262 


(45a) 


II 


The  inverses  are 


Z'' 


1.014975 
- .995025 


- .995025 
1.024876 


W = 


20.387249 

19.742660 

19.889737 


19.742660  19.889737 

20.830026  19.742660 

19.742660  20.387249 
_ 


20.436714  19.841465 

19.841465  20.239281 


The  rest  of  Equations  (13)  can  now  be  evaluated  as 


X = 


20.038895 

19.645015 

20.038895 


19.940670 

20.038892 

19.940670 


Y = X^. 


Equation  (16)  gives 


(45b) 

(45c) 

(45d) 

(46a) 

(46b) 

(47) 

(48) 


The  1st -order  residuals  are  calculated  as  in  Equation  (38): 


- .002 

- .026 


.043 

-.038 

.019 


(49) 


Their  standard  deviation  is 


a^(,)  = .0.33, 

which  is  but  slightly  larger  than  the  value  in  Equation  (42)  that  is  obtained  with  the  full  data  set. 


(50) 


It  is  interesting  to  see  the  effect  on  the  solution  of  “filling  in”  the  missing  data  point  in  Equation  (43). 
One  possible  procedure  is  to  replace  the  missing  point  by  one  half  the  sum  of  the  two  averages  on  the 
defective  row  and  column.  (For  a larger  data  grid  one  would  not  average  over  the  entire  row  and  column, 
but  only  over  points  near  the  missing  data  point).  Doing  this  we  get,  instead  of  Equation  (43), 
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a(0)  = 


1 6 

- 1.75  -2 
-4  1 


(51) 


Since  now  again  we  have  a full  data  set,  we  may  use  the  much  simpler  solution  of  chapter  3.  Keeping  the 
same  a-priori  variances  O/^  = = 10,  the  solution  is 


/ 3.458\ 

- 1.890  , 
\- 1.517/ 


J '-595^ 

1 .645 1 ’ 


(52) 


with  resulting  Ist-order  residuals 


a('»  = 


- .863 

.897 

1.736 

- 1.754 

-.888 

.872 

whose  standard  deviation  is 


a^(,)=  1.356. 


(53) 


(54) 


This  is  much  worse  than  the  corresponding  results,Equations  (49)and  (50),  where  the  missing  data  point 
was  not  filled  in  but  was  given  zero  weight  in  the  solution.  This  is  undoubtedly  dueto  the  very  small  data 
size.  For  larger  data  grids,  with  relatively  fewer  missing  data  points,  the  degradation  of  the  solution  due  to 
filling  in  the  missing  points  is  not  expected  to  be  so  severe. 


D.  EXAMPLE  4 

In  this  example  biases  as  well  as  bias  rates  are  considered.  To  each  of  the  biases  of  Example  1 we  add 
bias  rates  as  shown  in  Figure  2.  Remember  that  the  time  origin  t = 0 with  respect  to  which  the  bias  rate  and 
higher-order  effects  propagate  along  each  track  must  be  placed  at  the  middle  of  the  grid.  The  data  matrix 
computed  from  Figure  2 is 


A<0>  = 


I 5 
-6.5  -3 
-3.5  1.5 


which  has  the  standard  deviation 


Dl.  Bias  Solution 


°^(0)  -4.164. 


(55) 


(56) 


As  in  Example  2,  we  take  for  the  a-priori  bias  standard  deviations 

— 0^  — 1 0. 
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Figure  2.  Biases  and  Bias  Rates  for  F.xample  4 


The  solution.  Equations  (3 1 ),  is 


with  a standard  deviation 


V') 


.342. 


D2.  Bias  Rate  Solution 


The  Ist-order  residuals,  Equation  (58),  are  now  used  to  solve  for  the  bias  rates.  We  assume  a-priori 


values  on  bias  rates  of 


' = ( 


The  U and  V matrices  (see  Table  1)  are  readily  constructed  from  Figure  2. 


- .5 

- .5 

- -sl 

- 1 - 1 

u= 

_ .5 

.5 

■d- 

V = 

- .5  - .5 

1 1 

whence  F.quation  ( 1 2c)  yields 


Q = - 


1 

T 


2 

1 

.T 


(60) 


(61) 


Some  of  the  intermediate  results.  Equations  (12)  and  ( 1 3)  are : 


A=.54/3^3.  4/  = 2.29/2^2 

.321659 

-.109170  .218341 

= 

-.109170 

.485415  .109170 

■ 

.218341 

.109170  .321659 

Z"'  = 

'1.248333 

J. 04 1667 

1.041667"]  , 
1.248333J  ’ 

r 10.148762 

4.148438 

-8.296910" 

w = 

4.148438 

3.926063 

-4.148438 

-8.296910 

-4.148438 

10.148762 

z = 

2.637703 

-2.201022 

-2.201022 

2.637703_ 

■ 4.480301 

-4.480301" 

X = 

2.240150 

-2.240150 

-4.480301 

4.480301 

Y = XT. 


Equations  (16)  yield  the  vectors 


(62) 


(63a) 


(63b) 


(63c) 


(63d) 


(63e) 

(630 


(64) 
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whereby  the  solutions,  Equations  (14),  are  obtained  as  follows: 


Row  bias  rates: 


/ .299  \ 
b«)  = X = [-  .385  |; 


Column  bias  rates: 


b(C)  = V = 


.242) 


These  solutions  compare  very  badly  with  the  “true”  bias  rates  as  given  in  Figure  2: 


(65a) 


(65b) 


^TRUE  i 0 j ’ ^TOUE  (.2)  • 


(66) 


The  reason  is  that  the  amount  of  data  is  very  small.  With  only  two  intersections  per  row  and  three  per 
column,  one  would  not  expect  the  data  to  have  much  strength  for  bias  rate  determination.  Tfie  bias  solu- 
tion, Equation  (57)  took  out  too  much  of  the  residuals  in  this  case. 

Subtracting  the  effect  of  the  bias  rate  solutions  from  Equation  (58)  leaves  the  2nd-order  residuals 


.002  .025 

- .009  - .038 
-.012  .015 


which  have  a standard  deviation 


(67) 


(68) 


D3.  The  Correlation  Matrix 

It  is  instructive  to  look  at  the  correlation  matrix  for  this  example.  We  define  it  as 


(69) 


where  P is  the  usual  covariance  matrix 


= j?“l 


P = B 


(70) 


That  is,  the  diagonal  elements  of  R give  ihe  solution  variances,  while  the  off-diagonal  elements  are  the  cor- 
relation coefficients. 
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For  the  bias  solution.  Equation  (57),  the  correlation  matrix  is 


Rows 

Columns 

1 2 

3 

1 

2 

20.358  .976 

.976 

.984 

.984 

20.358 

.976 

.984 

.984 

20.358 

.984 

.984 

20.226 

.984 

20.266 

(71.) 


The  diagonal  elements  show  that  the  a-priori  uncertainty  of  o^  = 100  has  been  reduced  to  approximately  20. 
This  reduction  was  a little  more  effective  for  the  column  biases  than  for  row  biases,  which  is  to  be  expected, 
since  the  columns  have  more  intersection  data  than  the  rows.  The  correlation  coefficients  are  all  close  to 
unity  indicating  that  all  biases  may  be  in  error  by  the  same  amount. 


The  correlation  matrix  for  the  bias  rate  solution,  Equations  (65),  is  more  interesting: 


10.149  .657 

-.818  j 

.866 

-.866 

3.926 

-.657 

.696 

-.696 

10.149  1 

-.866 

.866 

1 

1 

1 

2.638 

-.834 

2.638 

Now  the  a-priori  variance  o^  = 25  has  been  reduced  to  approximately  10  for  the  first  and  third  rows,  to 
3.9  for  the  second  row  and  to  2.6  for  both  columns.  The  correlation  coefficients  are  not  close  to  unity 
anymore.  Thus,  the  element  /?  [3  = -.818  implies  that,  if  the  error  in  the  bias  rate  solution  for  the  first  track 
is  positive,  the  error  for  the  third  track  will  most  probably  be  negative.  Similarly,  all  the  other  elements  can 
be  qualitatively  understood  in  terms  of  the  geometry  of  Figure  2. 


E.  EXAMPLE  5 

The  previous  examples  had  only  five  intersecting  tracks.  They  are  intended  mainly  to  serve  as  check- 
out cases  for  a computer  program.  This  and  the  next  example  are  more  realistic.  They  consist  of  a total  of 
iOOO  intersections  with 

n = 50  rows,  m = 60  columns. 


The  data  matrix  was  constructed  by 


Af=h,-b,; 


1 < / < 50 
I < / < 60 


(7.1) 


where  the  110  biases  b,  and  bj  were  drawn  from  a Gaussian  distribution  with  zero  mean  and  standard 
deviation 

Of,  = 5.  (74) 


1 


A 
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Two  variations  of  this  data  grid  were  also  constructed.  The  first  had  random  noise  6,y  added  to  the 
This  noise  was  generated  from  a uniform  distribution  -.5  < fi,y  < .5,  so  that  its  standard  deviation  is 

Ofc  = -;!=  = .2W.  (75) 

The  second  variation  consisted  in  randomly  deleting  26%  of  the  data  in  liquation  (73). 

Instead  of  providing  the  a-priori  variances  and  as  input  values,  they  are  calculated  in  this  as  well 
as  the  next  example  by  a prescription  which  is  given  in  chapter  5. 

The  results  of  the  solution  for  this  example  are  given  in  Table  3 in  terms  of  the  standard  deviations 
a^{k)  for  the  Ath-order  residuals  A***.  The  a-priori  standard  deviations  and  a'/,,  and  the  solution  error 
standard  deviation  a^t,  are  also  given.  With  the  full  data  set,  and  no  noise,  just  the  Oth-order  solution  (biases) 
reduces  the  residual  standard  deviation  from  6.87  to  .0013.  The  solution  errors  Al>  have  a standard  devia- 
tion = .0(X)‘)5.  Additional,  higlicr-order  solutions  reduce  the  residuals  further  by  negligible  amounts. 
This,  of  course,  is  expected,  because  the  data  consist  of  biases  only.  The  third  column  of  Table  3 shows  the 
results  for  the  noisy  data.  The  solutions  arc  able  to  reduce  the  residual  variance  only  to  the  noise  level, 
Kquation  (75).  The  last  column  is  for  the  case  when  26%  of  the  data  have  been  randomly  deleted.  Here 
the  results  shove  only  a slight  degradation  as  compared  with  the  full  data  set. 


Table  3.  Standard  Deviation  a^(*)  of  Ath-Ordcr  Residuals  A**  > for  F.xample  5 


4- 

Full  Data. 

Full  Data 

26%  Data  Deleted, 

No  Noise 

_ _ . 

With  Noise 

No  Noise 

0 

6.86827 

6.87238 

6.86084 

1 

.00133 

.28229 

.00182 

■> 

.00133 

.27562 

.00182 

3 

.00133 

.27314 

.00181 

Solution  Error  r;^/, 

9.53*1  O''* 

.03769 

.00131 

A-Priori  Standard  I 
Deviations  | 

9.713 

1.22*  lO'** 
5.88*10'*’ 

9.719 

2.58*10'^ 

1.22*10'^ 

9.703 

1.66*10"* 

8.01*10'*’ 

F.  EXAMPLE  6 


This  is  the  same  as  the  previous  example  with  the  addition  of  bias  rates.  The  bias  rates  were  selected 
from  a Gaussian  distribution  with  zero  mean  and  standard  deviation  = 0.1 . The  data  grid  is  assumed  to  be 
uniform  with  unit  time  between  neighboring  intersections,  so  that  a bias  rate  = 0.1  yields  a total  change  of 
5 units  in  the  data  on  a column  and  6 units  on  a row.  The  affect  of  the  bias  rates  should  therefore  be 
approximately  equal  to  that  of  the  biases  (see  Equation  74). 


1 


1 


i 

i 

/; 

i 

a 


( 


. , i 

■i 

■ I 

Table  4 gives  the  results.  This  time,  as  expected,  it  takes  the  Ist-order  solution  (bias  rates)  to  reduce  | 

the  original  data  variance  to  an  acceptable  level.  The  full  data  set  with  noise  is  again  reduced  to  the  noise  | 

level.  Equation  (75),  but  the  deleted  data  set  now  shows  a much  larger  degradation  than  in  the  previous  | | 

example.  The  solution  errors  are  in  general  much  larger  than  in  the  previous  example  without  bias  rates, 
but  they  are  still  at  an  acceptable  level. 


Table  4.  Standard  Deviation  o^(*)  of  A:th-Order  Residuals  A***  for  Example  6 


it 

Full  Data,  No  Noise 

Full  Data  With  Noise 

26%  Data  Deleted,  No  Noise 

0 

7.27954 

7.28569 

7.26447 

1 

2.37727 

2.40096 

2.34391 

2 

.00379 

.27556 

.22535 

3 

.00379 

.27308 

.15971 

Solution  Error  ■ 

‘’a* 

.25521 

[^Ah 

.01 1 19 

A-Priori  Standard 
Deviations 

h 

10.295 

.218 

1.67- 10"^ 

10.304 

.220 

1.22-10-^ 

10.274 

.214 

9.94- lO'"* 

V.  A-PRlORl  VARIANCES 

The  a-priori  variances  af  which  are  needed  in  Equation  (10)  have  been  assumed  to  be  given  input 
parameters  in  most  of  the  previous  examples.  In  practice  this  will  be  the  case,  because  one  should  have 
some  fairly  accurate  a-priori  knowledge  about  the  uncertainties  of  the  orbit  biases  and  bias  rates.  If,  how- 
ever, such  knowledge  is  not  available,  the  method  given  here  can  be  used  to  estimate  the  a-priori  variances 
from  the  data. 

The  original,  Oth-order,  residuals  are  assumed  to  be  due  to  biases  only.  Equation  (2); 

Af=bi-b,:  (76) 

Assuming  these  biases  to  be  uncorrclated,  and  their  expectation  to  be  zero 

E(V/)  = o?5,^.  <77) 

E(6,.)  = 0,  (78) 

we  get  by  taking  the  expectation  of  the  square  of  Equation  (76); 

(79) 


The  first  order  residuals  are  modeled  by 


■ 


Again,  assuming  that  the  bias  rates  ft,  also  obey  Equations  (77)  and  (78),  we  have  upon  taking  the  expecta- 
tion of  the  square  of  this  equation; 


(81) 


where  r-  is  some  mean  over  rows  and  columns  of  the  squares  of  the  intersection  times. 
In  the  same  manner  we  derive  the  a-priori  variances  for  the  higher-order  solutions: 


(82) 


( 


f 

I 

f 

» 

\ 
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It  must  be  stressed  that  the  above  relations  are  recommended  for  use  only  when  there  is  no  independ- 
ent a-priori  knowledge  available.  In  particular  they  do  not  allow  different  a-priori  values  for  different  rows 
and  columns,  which  may  be  desirable  in  a practical  situation. 

The  a-priori  a’s  which  were  used  in  Examples  5 and  6 are  actually  twice  as  large  as  the  values  given  by 
these  e.xpressions.  We  found  by  solving  these  examples  with  different  a-priori  a’s  that  the  solutions  (biases, 
bias  rates  etc.)  are  not  sensitive  to  the  a-priori  values.  The  residuals,  however,  (excepting  the  noisy  data), 
are  quite  sensitive.  Doubling  the  a-priori  variances  gave  us  a little  better  residuals.  Another  reason  for  tak- 
ing sliglitly  larger  a-priori  variances  than  the  best  available  estimate  is  that  our  segmented  method,  because 
it  neglects  all  correlations  between  the  different-order  parameters,  will  always  produce  an  optimistic 
covariance  matrix  P,  Equation  (70).  Increasing  the  a-priori  variances  will  correct  this  error  somewhat. 

VI.  CONCLUSION 

A method  has  been  presented  for  extracting  biases,  bias  rates  and  higher-order  terms  from  intersecting 
satellite  altimetry  tracks.  This  method  consists  of  segmenting  the  least-squares  solution  such  that  at  each 
stage  only  one  set  of  terms  is  being  solved  for.  It  is  clear  that  an  error  is  committed  in  this  approach 
because  the  correlations  between  terms  of  different  order  are  not  taken  into  account. 

In  order  to  estimate  this  error.  Example  (>.  with  1 10  tracks  having  both  biases  and  bias  rates,  was 
solved  by  a minimum-norm  least-squares  method.^  In  that  solution  the  model  was  not  segmented;  i.e.,  the 
biases  and  bias  rates  were  solved  for  simultaneously.  That  method  guarantees  the  absolute  minimum 
residual  variance,  while  constraining  the  parameter  vector  to  have  the  smallest  length.  For  the  perfect  data 
case,  column  2 in  Table  4.  that  solution  gave  a residual  o = 9,3’  10"’^  (all  of  which  is  probably  due  to  round- 
ol'O  as  compared  to  our  .00379.  For  the  noisy  data,  column  3,  the  result  was  .27,‘'54  compared  to  our 
.27556.  Furthermore,  our  error  variances  and  o^j,  are  acceptably  small.  (In  the  particular  minimum- 
norm  solution  used  they  came  out  to  be  much  larger  than  ours).  This  shows  that  the  segmented  algorithm 
is  able  to  reduce  the  residuals  very  close  to  the  optimum  level.  We  have  also  shown  that  we  can  approach 
this  optimum  level  arbitrarily  close  by  increasing  the  a-priori  variances.  However,  in  order  to  obtain  realis- 
tic covariances,  it  is  recommended  to  keep  the  a-priori  values  consistent  with  actual  a-priori  knowledge. 


« 
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The  advantages  of  this  algorithm  far  outweigli  the  tlicoretical  objections  due  to  the  segmentation  of 
the  problem.  These  advantages  are  ( I ) tlexibility  of  the  algorithm  and  (2)  computer  cost  savings.  Flexi- 
bility is  obtained  because  one  can  increase  the  number  of  parameters  by  going  to  higher-order  solutions 
without  increasing  the  complexity  of  the  problem  and,  more  important,  without  increasing  computer 
storage  requirements.  As  an  example  ol  cost  .savings,  the  minimum-norm  solution  in  Fxaxnple  6 required 
1 25  sec  CPU  time  on  the  CDC  6700  computer,  whereas  the  segmented  algorithm  took  only  55  sec.  This 
ratio  will  continue  to  favor  the  segmented  algorithm  with  increasing  number  of  higher-order  parameters  as 
the  square  of  that  number. 

It  is  expected  that  at  the  end  of  the  GF.O.S-3  mission  there  will  be  available  from  5000  to  1 0000 
satellite  altimetry  tracks.  Clearly,  even  the  segmented  algorithm  would  not  be  able  to  handle  a job  of  this 
magnitude.  Rather  than  a global,  world-wide  solution,  we  therefore  recommend  to  solve  smaller  patches  at 
a time.  Continuity  of  the  biases  and  higlier-order  terms  between  patches  would  be  assured  by  requiring 
some  overlap  of  neigliboring  patches. 

An  example  of  the  use  of  this  algorithm  with  real  geodetic  data  can  be  found  in  Reference  4, 
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APPENDIX  A 


COMPUTER  FORMULATION  FOR  THE  LEAST-SQUARES  ALGORITHM 


1.  INTRODUCTION 

This  formulation  gives  step-by-step  instructions  for  coding  the  general  solution  (Equations  (14),  (12), 
(13)  and  (16))  in  component  form  rather  than  in  the  matrix  form  which  is  given  in  the  text.  Refer  to 
Figure  1 for  notation. 


2.  INPUT  DATA 


A.  n = number  of  rows 
m = number  of  columns 

= a-priori  for  rows 

= a-priori  for  columns 


k = Q:  biases 
A:  = 1 ; bias  rates 
k = 2-.  2nd-order 


geoid-height  discrepancy 
at  (;,  /)**'  intersection  • 

= XXX  (anything)  at  points  • 
where  no  intersection  exists 


1 </<n 
\<j  <>m 


B.  For  bias  rates  and  higher-order  solutions,  you  will  also  need  the  following  data: 
= /*  time  on  row 
7-(C)  _ j-th  gp  yth  column 

~ intersection  time  on  row 
7^*^*  = last  intersection  time  on  /'*  row 

= first  intersection  time  on /*  column 
= last  intersection  time  on  /**'  column 


Shift  these  times,  so  that  the  new  time  r s 0 is  at  the  middJe  of  each  track: 


= t(R) 

I 

/' 

n 

‘ 2 

1 

c. 

E-; 

II 

1 

i/ 

'/ 

2 

'IP 

II 

III 

:0 

at  all  intersections 
• where  A,-,-  exist 


3.  PRELIMINARY  CALCULATIONS 
The  Mean  of  Residuals: 


I / 


The  Residual  Variance: 


2-  1 
N-  1 


I I 


Note: 


a.  The  ranges  of  / and  / in  these  two  sums  are  only  over  those  points  • where  A,,  exist.  N = total  numbers 

of  such  points.  '' 

b.  The  same  formulas  will  be  used  for  the  means  A*,  and  variances  a|  for  the  higher-order  solutions. 

c.  Print  Aq  and  . 

4.  BIAS  SOLUTION 

A.  Construct  the  (m  X n)  matrix  U and  the  « X m matrix  V: 


Ujl  - I at  intersections  • ^ 

= 0 at  all  missing  points  ' 


V.^.  - 1 at  intersections 

= 0 at  * 


I < / < n 
I </<m 


A-2 


I 


r 


1 


Go  To  MAIN  ALGORITHM  =*  Xj,yi 
The  solutions  are: 

^(R)  _ Pjjj.  j-th  j ^ ^ „ 

= Vj  = bias  for /**’  column,  1 </  < w 


B.  Bias  Removal 

AyP  = - X,-  + y’j  at  all  intersections  •. 

Compute  and  print  A|  and  o^  (see  step  3). 


5.  BIAS  RATE  SOLUTION 


A.  Construct  the  matrices  U and  V : 

I < I < M 
1 < / < m 


{/..  = T^'> 
II  It 


V..  = 

a u 


Go  To  MAIN  ALGORITHM  =>  x„.(y 
The  solutions  are: 

ith 


b/  = Xj  = bias  rate  for  / row;  1 < / < m 

it/  - Tj  - bias  rate  for  /*  column ; 1 </  < m 
B.  Bias  Rate  Removal 


A^?^  = A^l^  intersections  • 


Compute  and  print  A2  and  02- 


6.  HIGHER-ORDER  SOLUTIONS, k>l 

A.  For  k'^' -order  (k  = 0,  biases,  and  k = 1,  bias  rates,  were  done  in  steps  4 and  5). 
Construct  the  matrices  U and  V: 
t/,.  = (r<«>)* 


1 < /■  < M 


y = ,.(CUk  f 1</<W 

a ii  ’ ' 


A-3 


I 


Go  To  MAIN  ALGORITHM  =» 

Xj  = it^-order  solution  for  i'*’  row;  1 < / < « 

= jt^-order  solution  for /**'  column;  1 </  < w 


B.  Remove  k^-order  effect 

^(*+1)  = ^(*)  .XjUjj  + VjVjj  at  all  intersections. 
Compute  and  print  and 


7.  THE  MAIN  ALGORITHM  (for  fc*  iteration) 


rn 


1 </<n 


";■£  ’'I 


O . =-  V.  U: 
1/  n 


1 <j<m 

1 </<M 

1 <i<m 


(H'-'),*=x,5,*-  E -;r.QijQkr 

(2'')//  = ^A7-E  yW,; 

,=  i ' ' / 

INVERT  the  above  symmetric  matrices  B'"’  and  Z”* 
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